Frequency-dependent transport through a quantum dot in the Kondo regime 
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We study the AC conductance and equilibrium current fluctuations of a Coulomb blockaded 
quantum dot in the Kondo regime. To this end we have developed an extension of the numeri- 
cal renormalization group suitable for the nonperturbative calculation of finite-frequency transport 
properties. We demonstrate that AC transport gives access to the many-body resonance in the equi- 
librium spectral density. It provides a new route to measuring this key signature of Kondo physics 
which so far has defied direct experimental observation. 
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In the past years semiconductor quantum dots have 
gained considerable attention as tunable magnetic impu- 
rities ll|. Due to their small size, electronic transport 
through these structures is strongly influenced by the 
Coulomb blockade |3| ■ Quantum dots with an odd num- 
ber of electrons display a well-known many-body phe- 
nomenon, the Kondo effect |3j |j], as predicted in [5|. 
In these systems, a single unpaired spin is screened at 
low temperatures, giving rise to an enhanced DC (zero- 
frequency) conductance at low bias voltage. 

Theoretical studies have so far mostly focused on DC 
transport 5*. From the differential conductance dl/dV 
at finite bias voltage V one can obtain information about 
the Kondo resonance in the spectral density. Such mea- 
surements on a typical, symmetrically coupled quantum 
dot, however, access only the spectral function under 
non-equilibrium conditions. It is known that in this case 
a finite bias splits and suppresses the Kondo peak 6] . It 
has also been demonstrated experimentally that a suffi- 
ciently strong AC modulation of the g^e or bias voltage 
likewise suppresses the Kondo peak [21 |3|- Thus, it is 
by no means straightforward to measure the equilibrium 
spectral function of a Kondo quantum dot. Previous pro- 
posals to do so require either very asymmetric couplings 
to the leads or quantum dots in a 3-terminal geometry 

ilia. 

In this Letter we point out that the Kondo resonance in 
the equilibrium spectral density of a generic quantum dot 
can be extracted from AC transport measurements. We 
derive a direct relation between the spectral function and 
the linear AC conductance through a quantum dot [Eq. 
(O below]. Alternatively, the Kondo peak can be mea- 
sured by detecting equilibrium fluctuations in the current 
through the quantum dot at frequencies of the order of 
the Kondo temperature Tk [Eq. (jH)) below]. Such mea- 
surements have come into experimental reach with recent 
advances in the measurement of high-frequency current 
fluctuations [Ij. Conductance measurements in the rel- 
evant frequency regime are challenging, mainly due to 
parasitic currents through hardly avoidable capacitances. 
One should, however, be able to extract the signal due 
to the current through the quantum dot from this back- 
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FIG. 1: Single-level Anderson model with time dependent 
chemical potential /ic,(i) = /i ± (el4c/2) coscijt in lead a £ 
{L,R}. 



ground by an adiabatic ch ang e (much slower than the in- 
verse Kondo temperature |lj| ) of gate voltage that drives 
the dot periodically into and out of the Kondo regime. 

We apply the nonperturbative numerical renormaliza- 
tion group (NRG) in combination with the Kubo formal- 
ism [ij, [ijj to predict the frequency-dependence of the 
AC conductance and the equilibrium noise. Our results 
are numerically essentiall y ex act, in contrast to previous 
perturbative calculations [15J , and generalize those of [ij] 
to finite frequency. This allows us to access the experi- 
mentally most relevant regime of frequencies of the order 
of the Kondo temperature for a device that exhibits a 
strong Kondo effect. 

The exact correspondence between spectral function 
and AC conductance holds only when charge fiuctuations 
on the quantum dot are negligible. Using our rigorous 
NRG results, we verify that this condition is well satisfied 
for frequencies below the charging energy of the dot. 

The setup that we study is shown in Fig. Q] At low 
energies, a Coulomb-blockaded quantum dot is well de- 
scribed by the single-level Anderson model 



H = ^€dadid„-i-Und^ndi + ^Va{dlc,^^^ + h.c.) 

+ Yl N-Aia(<)]cj'„^C^„<„ (1) 



^e{L,R} 



where the chemical potentials lJ'L{R)it) 



M ± 



(eVacf^.) COS ujt for the left (right) lead include a time- 
dependent bias Vac- We assume that the AC pertur- 
bation does not couple to electrons in the local level 
d^ {nj_„ — d\d^). Va couples the level d to electron 
states Ckaa with momentum k and spin ct = ±1 in lead 
a G {L, K). In the following, D denotes the conduction 
electron bandwidth, eda = Cd + uB* /2, B* = gfisB, is 
the energy of the local level, including the Zeeman shift 
in presence of a magnetic field B. U is the Coulomb 
repulsion of electrons on the dot. 

Following [13 we relate the current / through the 
quantum dot to the local Green function of the level d. 
For frequencies hco <C Ac much smaller than the charge 
excitation energy Ac = min{|ed|, |C/-|-ed|}, charge fluctu- 
ations on the quantum dot can be neglected and the cur- 
rents flowing through the right and the left lead are equal 
to a good approximation: I^ = —Ir- It is then advanta- 
geous to write the total current as / = (/l — A/_r,)/(1 + A), 
where A = Fi/Ffl, Fq, = tti^V^ and v is the conduction 
electron density of states. Expressing the currents II and 
Ir by the Green function of d (see Eq. (15) in [l£j) and 
Fourier transforming from time t to frequency a), we find 
for energy-independent couplings F^ the current 
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h Lu 1 + A ./_^ ./_^, ./ 2tt 



I de 
dh / TT- 



X Re e*'(*-*i)(sincjt-sinwti)/(e)g"(t,ii) (2) 

in linear response to T4c- Here /(e) is the 

equilibrium Fermi function of the leads and 
g''iti,t2) = -iO{ti-t2){{d{ti),d^t2)}) the retarded 
Green function of d evaluated in equilibrium (14c = 0). 
From Eq. Q we obtain the real (G") and imaginary 
(G") parts of the linear conductance. 
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(3) 
(4) 



where F = Fl + F^ and A{uj) = -ilme''(w) is the 
spectral function. At zero temperature kT <^ h]^] -C Ac 
Eq. Q is readily solved for the symmetrized spectral 
function ^(e) = [A{e) + A(-e)]/2, 



Asiuj) 



d 



e2 Fl Tr duj 



[uG'iuj)]. 



(5) 



In the Kondo regime {{nd) — 1), the Kondo resonance 
is centered almost exactly around 0, thus the Kondo 
peaks in As{ijj) and A{ijj) are essentially indistinguish- 
able. Hence, Eq. lO allows us to extract the Kondo peak 
in the equilibrium spectral function from a measurement 
of the AC conductance. 
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FIG. 2: (a) Real and imaginary parts of the AC conduc- 
tance obtained via the Kubo approach (solid lines) and from 
the conductance formulas Eqs. Q and (dashed Unes). 
Note the excellent agreement for frequencies Tiuj < 100 k Tk 
(Ac ~ GOOfcTK). (b) Doniach-Sunjic tails in A{uj) are re- 
sponsible for the power-law behavior (dashed line) in the con- 
ductance (solid line) which is well described by G' ex a;~^" 
for 2QkTK < w < 200 fcT^ ^. (c) Comparison of the 
frequency-dependent conductance phase (j}{ui) in the Kondo 
regime (solid line) (parameters U = 0.12D, U/T = 6, e^ = 
— C//2, T = 0) with that for a resonant level {td = 0) of width 
kTx (dotted hue). 



The frequency-dependent equilibrium current fluctua- 
tions (Johnson-Nyquist noise) 



Ciu;)^ 



dte'^' [(/(0)/(<)) - (/)2 



(6) 



are related to the linear conductance by the fluctuation 
dissipation theorem (FDT) 



C{co) 



2huj 



exp{huj / kT) — 1 



G'iu;). 



(7) 



Consequently, the spectral function can alternatively be 
inferred from a measurement of C{lu). At zero tempera- 
ture, kT <^ hliul <C Ac, we arrive at 
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(8) 



Note that at zero temperature C{uj) is non-vanishing 
only for a; < 0, that is fluctuations have to be mea- 
sured by probing absorption by the quantum dot. In 
|ll| the measurement of current fluctuations at frequen- 
cies up to tj ~ 100 GHz has been reported. This fre- 
quency scale is of the same order of magnitude as the 
Kondo temperature in typical Kondo quantum dots jl|. 
It makes noise measurements a promising candidate for 



experimental studies of the Kondo peak in the spectral 
function of a quantum dot. 

We turn now to a discussion of the frequency depen- 
dence of AC conductance and equilibrium noise in a 
Kondo quantum dot with typical experimental param- 
eters. We apply the numerical renormalization group 
(NR G), first used by Wilson to solve the Kondo prob- 
lem [13 ■ This method is nonperturbative and does not 
suffer from low-energy divergences common to scaling ap- 
proaches. In particular, it provides accurate results in 
the most interesting crossover regime hu ^ kT^- Within 
NRG we apply two independent approaches: Using the 
Kubo formula, we are able to calculate the conductance 
in the Anderson model Eq. (Q numerically exactly. On 
the other hand, using Eqs. Q and Q, we obtain conduc- 
tance and noise from the single-particle Green function 
of the d level, that we determine by NRG as well. By 
comparing the two approaches we shall demonstrate the 
validity of the approximation underlying Eqs. ^ and Q) 
for frequencies of the order of the Kondo scale. 

We apply the Kubo formalism following Izumida et 
al. |lJ|. We define an electric current from lead L to 

J „ 



R as I — 



(Nr) - {Nl) 



where iV„ = J2^^ c'^ 

is the total number of electrons in lead a and Na — 
jr[H,Na\. Introducing a linear response tensor a by 

{No) = aai3iuj)fi'f^, where Ml(_r) = ± (eT4c/2) coswt is 
the time-dependent bias applied to the left (right) lead, 
the total linear AC conductance takes the form 



G(cj) = — [aLLi^) + (Tjini^) - (tlr{^) - (^b.l{^)] ■ 

(9) 
The complex response tensor a can be written as 

<7ap{i~^) = i^ [Kapiuj) - Kap{Q)\, whcrc 



Kap{uj) 



^^^_St+^ut 



Nf,{0),N^it) ) (10) 



with S -^ 0+. (...) in Eq. H1(J|) refers to the equilib- 
rium expectation value with respect to the Hamiltonian 
H with T4c = 0. We evaluate this expression using 
NRG 13] . From the matrix elements of the current oper- 
ator {n\I\m) the imaginary part K'^^ is obtained directly 
while the real part K'^„ can be calculated via a Kramers- 
Kronig transformation. 

In Fig. [21 we compare the conductance obtained from 
Eq. (O with the result of the calculation in the Kubo 
formalism for e^ = —U/2 and B = Q. We find excel- 
lent agreement for frequencies below the charge gap Ac. 
Moreover the large frequency {huj ^ kTx) asymptotes 
of the conductance reveal a decay of G'{lu) ^ uj^^^^ 
[Fig. 121(b)]. This is expected as a consequence of the 
Doniach-Sunjic tails of the spectral function [l2| together 
with Eq. Q. The deviation of G'{uj) from the unitary 
limit (by ^ 6%) is due to systematic numerical errors 
accumulating in the NRG procedure. Fig. ^c) shows 
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FIG. 3: The equilibrium current fluctuations C(lj) depend 
linearly on uj for h juj < kTn 20], as shown for various mag- 
netic fields B. The inset shows that the fluctuations are spin- 
independent even in presence of a finite magnetic field B* 
{B* = 10 fcTif, (red) crosses: a —^, (black) solid line: a —i). 
Parameters: U = 0.12D, U/F = 6, ed = -U/2, T = 0. 



the frequency-dependent phase 0(w) of the linear con- 
ductance, Gitu) = lG(a;)|e*'^('^). 

The frequency-dependence of the equilibrium noise ob- 
tained by a direct numerical evaluation of Eq. ® within 
NRG is plotted in Fig.O As one expects, the fluctuations 
reach a maximum for frequencies huj ^ U. The linear 
behavior C{u}) ~ w at small w is a manifestation of the 
Fermi-liquid nature of a screened Kondo impurity. Fig. (31 
also shows how a finite magnetic field B suppresses low- 
frequency equilibrium fluctuations. The inset demon- 
strates that at low frequencies the noise for spin up and 
spin down electrons is identical, even in the presence of 
a magnetic field. This is because an impurity with large 
charging energy Ac :^ F is half-filled: (n^-f ) + {udi ) = 1 • 
As a consequence of the Friedel sum rule |l8| , the spec- 
tral densities at the Fermi energy ep for spin-up and spin- 
down electrons Ay^{eF) oc sin {TT{ndyi)) are then equal. 

We now use Eq. ((SJ to extract the spectral function 
As from the numerical current noise data of Fig. where 
charge fluctuations on the dot have been taken into ac- 
count. The results are shown in Fig. ^ They nicely 
demonstrate the splitting and the suppression of the 
Kondo peak upon increasing the external field B. A com- 
parison with the spectral function directly calculated by 
NRG confirms that Eq. (jSJ does indeed work very well 
for frequencies fnu <^ Ac. In particular, the Kondo peak 
in the symmetrized spectral function As{lo) obtained by 
our method is essentially indistinguishable from that in 
the spectral function A{uj) directly calculated by NRG. 
Due to intrinsic broadening in the NRG which affects 
both A{u!) and C(w), Eq. © becomes less accurate for 
finite B. 

Having demonstrated that neglecting charge fluctua- 
tions in deriving Eqs. © and ^ is very well justified. 
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FIG. 4: Spectral function A{uj) = As{u)) as extracted from 
the current noise C{lo) via Eq. Q (lines) compared to A{uj) 
obtained directly by NRG (stars) for different values of the 
magnetic field B. Eq. Q captures the Kondo peak in A{ui) 
very well, including suppression and splitting in presence of a 
magnetic field B. 



we comment on another limitation of our approach. Cal- 
culating the linear response of the quantum dot, we have 
assumed the bias voltage in a conductance measurement 
to be small. While this in itself is not a problem, one 
might worry, however, that due to the finite frequency of 
Vac extra decoherence processes would make the limit of 
linear response very restrictive. The decoherence rate r 
of the dot's spin due to the oscillating bias voltage for 
huj ^ k Tk can be estimated |lj| as 



TiT- 



eK, 



kT, 



K 



kTK \kTKj fiuj [\n{hujlkTK)Y 



(11) 



For the Kondo physics not to be disrupted by T4c, we 
need h,T~^ /kTK «C 1. Eq. l(TT|l shows that this con- 
dition can be easily fulfilled by the usual requirement 
eVac ^ kTx- We expect that this statement remains 
true for frequencies of the order of Tk- We conclude 
that transport is described accurately by our linear re- 
sponse conductance for voltages that are much smaller 
than the Kondo temperature. 

In conclusion, we have studied AC transport through 
a quantum dot in the Kondo regime, using the numeri- 
cal renormalization group technique in combination with 
the Kubo formalism. We have expressed linear conduc- 
tance and equilibrium current fluctuations in terms of the 
single-particle Green function. This relation becomes ex- 
act at low frequencies, when charge fluctuations on the 
dot can be neglected. It has been shown to work very 
well for frequencies of the order of the Kondo scale. This 
opens up the exciting possibility of measuring the equilib- 
rium Kondo resonance directly in a transport measure- 
ment. 
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